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Abstract: Recursive solution for a general homogeneous in time pure birth branch- 
ing process with simultaneous production of any number of "particles" and with 
continuous evolution parameter is given. Calculational algorithm based on the use 
of Koenigs function and functional Schroder equation is described. It is shown that 



multiplicity distributions in e + e~ annihilation into hadrons for cm. energies up to 
189 GeV are well described by the modified negative binomial distribution, explained 
by simple pure birth branching process without multiple simultaneous "particle" pro- 
duction. The energy dependence of the evolution parameter is also discussed. 



1. Introduction 



Recently it has been shown [1-13] that multiplicity distributions in e + e~ annihilation 
into hadrons are fairly well described by the modified negative binomial distribution 
(MNBD). Two scenarios based on the use of the pure birth branching processes 
(for description of the mathematical formalism see for example fl4| , [15]) have been 
proposed to explain the origin of the MNBD. In the scenario advocated in the pa- 
per [H the binomial distribution of particle production sources occurs at some initial 
stage, each source evolves afterwards according to the pure birth branching process 
with immigration. In other scenario || || a fixed number of particle production 
sources occurs initially and they evolve also according to the pure birth branching 
process with production of intermediate neutral clusters. The final state hadrons 
are produced through decay of these clusters. In the branching process used in both 
scenarios no multiple simultaneous "particle" or cluster production is allowed. In 
this letter we study the possibility of having simultaneous multiple "particle" pro- 
duction. In section 2 a recursive solution for the pure birth branching process with 
simultaneous production of any number of "particles" is given together with a cal- 
culational algorithm, based on the use of Koenigs function and functional Schroder 
equation|l6| (recent description can be found in [|17], [18|]). In section 3 results of fits 
of the charged particle multiplicity distributions in e + e~ annihilation into hadrons at 
cm. energies above 20 GeV [19-28] to the MNBD and to the distribution resulting 
from the branching process with additional simultaneous pair "particle" production 
are presented. Our conclusions and discussion are given in last section. 



2. Recursive solution for a general pure birth branching pro- 
cess 

Let us recall that the branching process with continuous evolution parameter t is 
determined by the differential probability densities ai for the transition of one ob- 
ject, "particle", into i objects, "particles". All "particles" are assumed to evolve 
independently. For a homogeneous branching process the oti do not depend on t and 
probabilities for transition 1 — > i in an infinitely small interval At are 

p 1 ^ i (t,t + At) = aiAt ,2 = 0,2,3,... (2.1) 
and for transition 1 to 1 



p 1 ^ 1 (t,t + At) = 1 - a At 



[2.2) 



with 



;2.3) 



1 



For the pure birth branching process ao = 0. 

The probability distribution Pi(t) = pi^j(0, t) for the process having one particle 
at t — can be found from the forward Kolmogorov equation 

din Otti 

— — = fix)—— with initial condition Tn\t=o — % , (2-4) 
at ax 

for the probability generating function 

oo 

m{x,t) =J2Pi(t)x i , (2-5) 
i=i 

with 



f( x ) = aiX% ~ ax ■ (2-6) 

i=2 

Using the Taylor expansion of the equation [O] over x one can obtain the following 
system of differential equations for the probabilities 



dp 

dp 2 



dp 3 
dt 



(jf Oi 2 p\ - 2ap 2 , (2.8) 

a 3 pi + 2a 2 p 2 - 3ap 3 (2.9) 

(2.10) 



and for arbitrary k 



dpk 
dt 



fc-i 

joik-j+iPj - kap k (2.11) 

3=1 



with initial condition 



Pfc(0) = 5 lk . (2.12) 
These equations can be solved one after another, for p\ 

Pi = exp (—at) , (2-13) 

for p 2 

p 2 = a 2Pl (l - Pl ) , (2.14) 

where ai = ai/a , and so on. 

One can observe that pk has polynomial dependence of fe-th power on p\ 



2 



(2.15) 



i=l 



with constant coefficients Kik- These coefficients can be found using the following 
recursions 



TTll = 1 , 7Ti2 = -7T 2 2 = Gt 2 , 



(2.16) 



fc-1 fc-1 
(k-i)7Tik = X^ a (fc+i-i) 7I X? = ZK^ _ ^') a (i+i) 7r ^-i) = 1,2,..., A; - 1 , (2.17) 



and the coefficient Tr^ can be found from the initial condition 2.12 



fc-i 



(2.18) 



1=1 



This solution satisfies also the backward Kolmogorov equation dm/dt = f(m). 
Solving the backward equation in quadratures one can obtain a functional relation 
between m and x. For example, for the branching process with non-zero a 2 and 
one has: 

(p(m) = <p(x) exp(-at) (2-19) 



with 



and 



(p(x) 
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x 



. . i , . i 

(1 — x)^+ l [x + 7)^ 

a 2 + a 3 a 



(2.20) 



(2.21^ 



«3 



«3 



In practical calculations the polynomial formulae |2.15| for pk have limited nu- 
merical precision. More robust calculational algorithm based on the mathematical 
formalism used by Fatou [|i6| in the studies of the iterations of rational functions is 
desribed briefly below. Let us denote 



m n (x,t) = m(x,nt) 



(2.22) 



The following relations are valid for the probability generating functions m n (x,t) 
(see for example |L4|, PH ) 



m n (x,t) = m n -i(m(x,t),t) = m(m n -i(x, t), t) 



(2.23) 



In analogy with |16j one can introduce the Koenigs function for the iterations |2.23| 
as a limit: 
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K(x) = lim 



m n (x,t) 



(2.24) 



For the pure birth branching process the Koenigs function has the following form 



(2.25) 



8=1 



Let us denote Tin = K{. The recurrence relations |2 . 1 7| lead to the following recurrence 
for Ki 



i-1 



i-i 

-3 = yi-l"'' J - '- ''J 

i=i 



(2.26) 



The remarkable property of the Koenigs function [16], [L7|, [Ll| is the functional Schroder 
equation: 



K (m(x, t)) — p\K(x) , 
it follows as a limit from the functional relations |2.23 L . Let us denote 



rrv 



x J 



,J = 2,3, 



(2.27) 



(2.28) 



The use of the Taylor expansion of the Schroder equation |2.27| results in the following 
relations: 



U) 



i = 1,2,3,... 



i=i 



And can be found using successively the convolutions: 

i-i 



Pf = J =2,3,... 



vT ) = Y,W 1 j-x j = N,N + l,... iV = 3,4, 



i=i 



(2.29) 



(2.30) 



(2.31) 



Using [2.29| , |2.30| and |2.31| one can obtain the recurrence for p n . Indeed, let us assume 
that we know Pi for i — 1,2, . . . ,n — 1. Then at step n we calculate 



n— j — 1 

p ti) - V" „.„'.'- u 



£ PiPn-i J — 2, . . . , n — 1 



(2.32) 



i=i 



1 It is worth to note that the Schroder equation and the equation 2.19 are identical, therefore the 



Koenigs function for the branching process with non-zero oti and coincides with <p(x) defined in 



2.2C. In general case the K(x) is proportional to exp(— a J y^) 
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and 



Pn — Pi 



(2.33) 



Finally, the equation 2.29 leads to the following expression 



n-l 



Pn = «n(Pl ~Pl)-J2 K iP 
J'=2 



Cj) 

n 



The probability to produce n_ negatively charged particles (n. 
scenario [||, || can be found using the formula: 



(2.34) 
n c h/ '2) in the 



oo oo 



Here iV denotes the number of initial sources of particle production, is the usual 
binomial probability to have n_ negatively charged particles ( we remind that n_ in 
this scenario is equal to the number of produced charged particle pairs) from decay 
of j clusters and e is the probability of cluster decay into pair of charged particles. 
In practical calculations we stop summation in |2.35| when the difference between 1 



and the sum of Pj becomes smaller than the computer precision. 
3. Results of fits 

The pure birth branching process with non-zero «2 and 03 is determined by these 
differential probability densities and by the evolution parameter t. The probability 
distribution pi(t) depends on some combinations of the a 2 and a 3 and on pi(t) (|2.13| ). 
The final multiplicity distribution (|2.35|) is characterized also by the number of initial 



sources N and by the probability e. As in the papers || | we fix iV at the value 
7. We use as free parameters in the fits pi, e and the ratio a = 0:2/(02 + 03)- The 
parameter a is equal to one for the MNBD and is equal to zero for the branching 
process with simultaneous pair "particle" production. Results of the fits are given 
in Table 1. In the cases when the three-parameter fit has minimum with a — 1, only 
the results of the two-parameter fit with fixed a = 1 are shown. 

The influence of the parameter a is illustrated in Figure [TJ, where the multiplic- 
ity distribution at y/s = 189 GeV is compared with the predictions of the three- 
parameter fit and with the predictions of the two-parameter fit with fixed a = 0.8, 
having p x = 0.509 ± 0.007, e = 0.861 ± 0.012 and x 2 /NDF = 19.1/(29 - 2). From 
Figure [I] it is seen that the predictions of the fit with a = 0.8 exceed the predictions 
of the three-parameter fit at n ch below 12, the opposite is observed at n ch above 40. 

One can see from Table 1 that the quality of the fits is quite good and that the 
values of e and x 2 f° r the MNBD fits are practically the same as in the previous 



5 



1 



i 



Vs = 1 89 Ge^ 



- 



2 



O data rj 
A fit with a = 0.987 ( J 

■ fit with a = 0.8 



i 



,5 



o 



20 30 40 50 60 



Figure 1: Multiplicity distribution measured by the OPAL Collaboration at yOj = 
189 GeV compared with the predictions of the fits with a = 0.987 and a = 0.8. 



publications [0, |[ (5|. The values of the parameter a are concentrated near one and 
the allowed fraction for simultaneous pair production is below a few per cent. We 
have checked also that the parameter a is near one in the fits with N different from 
seven, in most cases these fits have higher x 2 than the fits with N = 7 (see also 0). 
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exp. 


V 


r 1 


a 


e 


x 2 


N„ 


TASSO [19] 


22 


0.773±0.008 


1 


0.626±0.007 


22.1 


14 






0.854±0.025 


0.706±0.138 


0.658±0.138 


11.4 




HRS [20] 


29 


0.758±0.006 


1 


0.697±0.006 


7.8 


14 


TASSO [19] 


34.8 


0.688±0.004 


1 


0.668±0.004 


15.5 


18 






0.714±0.011 


0.948±0.024 


0.681±0.007 


9.8 




TASSO [19] 


43.6 


0.645±0.006 


1 


0.693±0.007 


27.1 


19 


AMY [21] 


50 


0.658±0.018 


1 


0.766±0.017 


2.3 


19 






0.674±0.068 


0.973±0.112 


0.772±0.029 


2.2 






52 


0.634±0.014 


1 


0.753±0.015 


5.8 


19 




55 


0.637±0.015 


1 


0.767±0.015 


4.6 


19 






0.650±0.042 


0.975±0.081 


0.772±0.022 


4.5 






56 


0.623±0.012 


1 


0.757±0.014 


11.0 


19 




57 


0.614±0.014 


1 


0.767±0.014 


7.5 


19 






0.616±0.038 


0.994±0.070 


0.768±0.021 


7.5 






60 


0.595±0.016 


1 


0.749±0.017 


5.8 


20 




60.8 


0.644±0.012 


1 


0.791±0.012 


14.4 


20 




61.4 


0.657±0.036 


0.935±0.075 


0.793±0.018 


9.4 


20 


ALEPH [22] 


91 


0.507±0.011 


1 


0.757±0.017 


7.4 


26 


DELPHI [23] 


91 


0.512±0.003 


1 


0.774±0.004 


30.1 


25 


L3 [24] 

L J 


91 


0.512±0.010 


1 


0.751±0.014 


11.8 


23 


OPAL [25] 


91 


0.501±0.005 


1 


0.776±0.008 


5.0 


25 


OPAL [26] 


133 


0.456±0.019 


1 


0.762±0.031 


5.3 


25 






0.525±0.059 


0.889±0.106 


0.821±0.054 


3.6 




OPAL [27] 


161 


0.441±0.017 


1 


0.768±0.031 


3.6 


25 


OPAL [28] 


172 


0.463±0.019 


1 


0.838±0.034 


5.8 


29 


OPAL [28] 


183 


0.437±0.012 


1 


0.826±0.020 


13.8 


29 






0.496±0.034 


0.904±0.060 


0.882±0.033 


10.3 




OPAL [28] 


189 


0.402±0.008 


1 


0.772±0.014 


3.3 


29 






0.411±0.025 


0.987±0.037 


0.780±0.026 


3.1 





Table 1: Results of the fits of the distribution 2.35 to the charged particle multiplicity 
distributions in e + e~ annihilation into hadrons. 



4. Discussion and conclusions 

In recent publication |L3[ Biyajima and coworkers have observed that the LEP data 
at -y/i = 133 GeV |26| favour logarithmic dependence ~ log y/s of the evolution 
parameter t in contradiction with log log Q 2 prediction given by QCD. For logarithmic 
dependence the probability pi should have power law dependence on the cm. energy 
y/s. The energy dependence of the parameter pi for the MNBD fits is shown in 
Figure^. From Figure it is seen that recent LEP measurements at 161 |27| 172, 183 
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Figure 2: Energy dependence of the probability p\ for e + e annihilation into hadrons. 
Solid line shows the predictions of the power law fit. 



and 189 p8| GeV also favour the logarithmic dependence of the evolution parameter 
on y/s. The data are fairly well described by the power law dependence pi = A(yfs)~ B 
with the slope B = 0.310 ± 0.005 and with x 2 /NDF = 104.1/(21 - 2). As noted 
in ||13||, this contradiction with QCD is probably attributed to hadronization effects. 



Our conclusions are the following. The recursive solution for the general pure 
birth branching process is given. The calculational algorithm based on the use of 
the Koenigs function and the functional Schroder equation is described. The results 
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of the fits to the charged particle multiplicity distributions in e + e~ annihilation 
into hadrons for cm. energies up to 189 GeV show the validity of the MNBD 
parametrization and the absence of the component with simultaneous pair production 
at the level exceeding a few per cent. The latest LEP measurements favour the 
logarithmic dependence of the evolution parameter t as noted earlier flT3|. 
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